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EXPRESS MAIL NO. EV335518195US 
METHOD AND SYSTEM FOR WATER QUALITY TRADING 

TECHNICAL FIELD 

[0001] The described technology relates generally to analysis of pollutant loading 

in a watershed and particularly accounting for the uncertainty in both the pollutant 
loading and pollutant removal efficiency of best management practices 
implemented as controls and in determining trading ratios. 

BACKGROUND 

[0002] Many watersheds throughout the country are approaching or already 

exceeding their allowable pollutant loading capacities referred to as total 
maximum daily loads ("TMDLs"). As a result, it is becoming increasingly 
important for watershed management to provide an economical way for ensuring 
that the TMDLs are not exceeded. Watershed management is important for 
watersheds not only because TMDLs for many watersheds are currently defined, 
but also because many existing TMDLs are expected to be revised and TMDLs for 
many other watersheds are expected to be promulgated in the near future. The 
changes to the TMDLs will likely have a significant impact on the treatment costs 
for both point and non-point sources of pollutants. The United States 
Environmental Protection Agency published a Water Quality Trading Policy on 
January 13, 2003, which Is hereby incorporated by reference, stating that a 
market-based approach to water quality trading can provide flexibility and achieve 
water quality and environmental benefits greater than would otherwise be 
achieved under more traditional regulatory approaches. A market-based, water 
quality trading approach can be used by utilities, industries, county governments, 
and so on to achieve economically sound development and meet the EPA's water 
quality goals. 
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[0003] Any such water quality trading should factor in the high degree of 

uncertainty that results from various watershed processes such as rainfall runoff, 
soil erosion, sediment delivery, pollutant generation, and so on. It is desirable to 
have a tool to assist in analyzing water quality trading potentials that factor in this 
uncertainty. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0004] Figure 1 illustrates an initial display page of the WQTS in one embodiment. 

[0005] Figure 2 illustrates a sample report listing the current and future pollutant 

loads for a selected confluence point in one embodiment. 
[0006] Figure 3 illustrates a display page for specifying the BMPs and treatments 

in one embodiment. 

[0007] Figure 4 illustrates a report listing the local credits for a subwatershed in 

one embodiment. 

[0008] Figure 5A illustrates a graph of the earned credits for the subwatersheds 

contributing to a confluence point in one embodiment. 
[0009] Figure 58 illustrates a report of the earned credits for a subwatershed 

contributing to a confluence point in one embodiment. 
[0010] Figure 6 illustrates a display page of a GIS interface for the WQTS in one 

embodiment. 

[0011] Figure 7 illustrates a map overlaid with the loading information of the 

subwatersheds in one embodiment. 
[0012] Figure 8 illustrates a map overlaid with the local credit information of the 

subwatersheds in one embodiment. 
[0013] Figure 9 illustrates a map overlaid with the earned credit information of the 

subwatersheds in one embodiment. 
[0014] Figure 10 is a block diagram illustrating the components of the WQTS in 

one embodiment. 

[0015] Figure 11 is a flow diagram illustrating the processing of a component that . 

identifies subwatersheds in one embodiment. 
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Figure 12 is a flow diagram illustrating the processing of a component that 
calculates the loading for the subwatersheds that are upstream from a confluence 
point in one embodiment. 

Figure 13 is a flow diagram posted in the processing of a component that 
calculates the local credits for the subwatersheds that are upstream from a 
confluence point in one embodiment. 

Figure 14 is a flow diagram illustrating the processing of a component that 
calculates the earned credits for subwatersheds that are upstream from a 
confluence point in one embodiment. 

DETAILED DESCRIPTION 

[0019] A method and system for calculating trading ratios between a buyer and a 

seller of water quality pollutant load credits is provided. In one embodiment, the 
water quality trading system ("WQTS") inputs a definition of watershed and land 
uses of the watershed. The WQTS then calculates the load credits (e.g., the 
difference between the total maximum daily load ('TMDL") of the pollutant and the 
actual load) of applying various controls (e.g., best management practices 
("BMPs") and treatments) to the point sources and non-point sources of the 
watershed. The load credits take into consideration the uncertainty of the 
calculated loads. The uncertainty can result from uncertainty in the export 
coefficient of a land use, in the removal efficiency of a BMP, and so on. The 
WQTS calculates an uncertainty ratio that will reduce the number of transferable 
credits generated by calculated non-point source reduction to account for 
variability in the effectiveness of the controlling practice (e.g., BMPs) 
implemented. For example, a seller may have a daily load credit of 500 lbs/day 
and the uncertainty is 0.7, then the available credits for sale are 500 x 0.7 = 350 
lbs/day. To account for the difference between the uncertainty ratios of the 
available credits of a seller and the needed credits of a buyer, the WQTS 
calculates a trading ratio to indicate the equivalency between the available and 
needed credits. In one embodiment, the trading ratio is the exchange ratio of the 
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[0017] 



[0018] 



seller divided by the exchange ratio of the buyer. When the exchange ratio of 
both the buyer and the seller are the same (e.g., both .7), then the trading ratio is 
1 , which means the available credits can be traded on a one-for-one basis with 
the needed credits. This happens mostly when the sellers and buyers are located 
in the same vicinity. When the trading ratio is greater than one, then an available 
credit is worth more than a needed credit. For example, if the trading ratio is 2, 
then the buyer would need to purchase 2 available credits from the seller for 
every needed credit. In contrast, when the trading ratio is less than one, then an 
available credit is worth less than a needed credit. For example, if the trading 
ratio is 0.5, then the buyer would need to purchase 1 seller credit for every 2 
needed credits. The WQTS thus provides a way to calculate the equivalency 
between available credits and needed credits so that potential credit trades can 
be evaluated. 

[0020] The WQTS allows a user to evaluate the feasibility of conducting pollutant 

trades while accounting for the economics of the trade. The WQTS models the 
pollutants of a watershed taking into account the contributing pollutants and their 
loads, the contributing sources and their locations, the uncertainty in the pollutant 
loads, and BMPs for non-point sources and associated uncertainty. The WQTS 
calculates credits and cost differentials among various sources and estimates 
trading ratios for buyers and sellers within a watershed. 

[0021] The WQTS initially calculates the pollutant loading of the subwatersheds of 

a designated watershed for use when analyzing the effects of various controls. In 
one embodiment, the WQTS uses a database that contains a description of river 
basins that specifies each confluence point of a river basin and each watershed 
and each subwatershed. After a user selects a confluence point at which the 
trading analysis is to be performed, the WQTS identifies upstream subwatersheds 
that contribute to that confluence point. The WQTS then calculates the pollutant 
load of each contributing subwatershed so that the effect of the controls can be 
analyzed. The database includes for each subwatershed a description of each 
current and planned land use (including the area and export coefficient of the 
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pollutant) and the load of each current and planned point source. In one 
embodiment, the WQTS calculates the total current and future pollutant loads for 
the point and non-point sources of a subwatershed. The current and future 
pollutant loads for point sources take into consideration the current point sources 
and the planned point sources (e.g., add a water treatment plant or remove the 
secondary treatment of an existing water treatment plant) for a subwatershed. 
The current and future pollutant loads for non-point sources take into 
consideration the current and planned land uses in the subwatershed. The 
WQTS also calculates an estimated future load at 50% and 95% confidence 
levels. 

[0022] The WQTS allows a user to specify the modifications (e.g., adding, 

removing, or changing) of the controls such as the BMPs of the land uses and the 
treatment of the point sources of a subwatershed. The WQTS allows a user to 
specify the land uses of a subwatershed and the percent of the subwatershed 
affected and to specify the BMPs that are to be applied to each land use, the year 
of installation, and whether the BMPs are to be in serial or in parallel. The WQTS 
also allows a user to specify the types of treatments along with the year of 
installation that are to be applied to point sources of a subwatershed. 

[0023] After the BMPs and treatments are defined, the WQTS calculates the local 

load credit of each subwatershed. The WQTS calculates the load of each land 
use of a watershed taking into consideration the percent of the subwatershed 
affected by the land use, the area of the land use, the export coefficient of the 
land use, and the pollutant removal efficiency of each BMP. The total load for the 
subwatershed is the sum of the load of its land uses and its point sources. The 
"local credit" of a subwatershed is the difference between the load prior to and 
after application of the controls and represents the cumulative effect of the 
controls of a subwatershed. 

[0024] When credits are to be traded, the WQTS calculates the "earned credit" of 

a subwatershed that represents the portion of the local credits that are realized at 
the confluence point taking into consideration the uncertainty. Various factors 
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influence what percent of a local load contributes to the load at the confluence 
point. The aggregate of the factors is referred to as the "exchange ratio." The 
load of a subwatershed that contributes to a confluence point is the exchange 
ratio times the local load, and the earned credit of a subwatershed at a confluence 
point is the difference between the load prior to and after application of the 
controls multiplied by the exchange ratio. The exchange ratio may aggregate the 
effects of the "uncertainty ratio," "delivery ratio," and so on. The uncertainty ratio 
represents the uncertainty in the calculated load of a subwatershed resulting from 
the variation in various parameters including the export coefficient of a land use 
and the removal efficiency of a BMP. The uncertainty ratio represents the percent 
of the load that will be generated with a certain confidence level. For example, if 
the load is calculated to be 100, then the load with a 50% confidence level may be 
90 and the load with a 95% confidence level may be 75. The delivery ratio 
represents the fraction of the load of a subwatershed that is actually delivered to 
the confluence point. Each downstream subwatershed on the way to the 
confluence point may only deliver a certain percent of its inflow to its outflow, 
which is represented by its delivery ratio. Thus, the effective delivery ratio of a 
subwatershed to a confluence point is calculated based on the delivery ratios of 
each downstream subwatershed. As discussed above, the WQTS calculates a 
trading ratio based on the exchange ratio for the available earned credits of a 
seller and the needed earned credits of a buyer at a certain confidence level. The 
exchange ratio can be used to determine the relative value of the available and 
needed credits. 

[0025] In one embodiment, the WQTS provides the cost of various BMPs for land 

uses and treatments for point sources. A user can specify various combinations 
of BMPs and treatments and analyze their cost-effectiveness. The WQTS in one 
embodiment may automatically identify an optimum combination of BMPs and 
treatments based on various constraints specified by a user. The optimum 
combination can be identified using various well-known techniques such as linear 
programming. 
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[0026] In one embodiment, the WQTS provides a geographic information system 

("GIS") that allows a user to identify a river basin and a confluence point of the 
river basin that is to be analyzed. The GIS database describes the river basins, 
the confluence points, the subwatersheds, the land uses (current and future), and 
the point sources. A river basin table provides geographic information of the river 
basin. The confluence point table provides for each river basin the location of 
each confluence point along with the "from" subwatersheds whose outflows are to 
the confluence point and the "to" subwatersheds whose inflows are from the 
confluence point. The subwatershed table provides geographic information for 
each subwatershed. The current land use table and the future land use table 
describe each current and future land use of a subwatershed and include its area 
and export coefficient. The point source table provides the current and future 
loads of each point source. 

[0027] The WQTS includes other database tables that describe the BMPs, the 

local delivery ratios, the confidence levels, the treatment types, the export 
coefficient and removal efficiency variations, and the TMDL at each confluence 
point. The BMP table provides the removal efficiency, the life, the cost, and the 
coefficient of variation of the removal efficiency for each BMP. The delivery ratio 
table provides the delivery ratio for each subwatershed. The confidence level 
table provides the 2 value of the standard normal distribution at various 
confidence levels. One skilled in the art will appreciate that different types of 
distribution can be used. The treatment type table provides the concentration of 
the pollutant, the life, and the cost of various treatment types. The export 
coefficient variations table provides the coefficient of variation of the export 
coefficient of each land use. The TMDL table provides the TMDL of the pollutant 
at each confluence point. 

[0028] Figure 1 illustrates an initial display page of the WQTS in one embodiment. 

The display page 100 includes a river basin field 101, a confluence point field 
102, and a calculation level field 103. The river basin field provides a drop-down 
list of the river basins defined in the database and allows a user to select the river 
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basin to be analyzed. The confluence point field provides a drop-down list of the 
confluence points within the selected river basin as defined by the database. The 
user selects the confluence point at which the pollutant load is to be analyzed. 
The calculation level field specifies the level (e.g., watershed or subwatershed) at 
which the analysis is to be performed. (This description assumes that the user 
has selected the subwatershed calculation level.) The display page also includes 
a loading button 104, a controls button 105. a trading button 106, and a 
requirements button 107. 

[0029] After a user selects the river basin, the confluence point, and the 

calculation level, the user selects the loading button to direct the WQTS to 
calculate the current and future pollutant loads for subwatersheds that contribute 
to the selected confluence point. The WQTS generates a report that lists for each 
subwatershed the total current load and total future load for the point and non- 
point sources. Figure 2 illustrates a sample report listing the current and future 
pollutant loads for a selected confluence point in one embodiment. The report 
also lists the total future loads with a 50% and 95% confidence level for non-point 
sources. The report lists the loads for each subwatershed and then provides the 
totals at the confluence points. 

[0030] After a user selects the river basin, the confluence point, the calculation 

level, and then the loading button, the user selects the controls button to specify 
the BMP to be applied to the land uses of a subwatershed and the treatments to 
be applied to point sources of the subwatershed. Figure 3 illustrates a display 
page for specifying the BMPs and treatments in one embodiment. The display 
page 300 includes a subwatershed field 301, non-point source area 310, and 
point source area 320. The subwatershed field identifies the subwatershed 
currently selected. Scroll buttons 330 allow a user to select different 
subwatersheds that are upstream from the selected confluence point. The non- 
point source area includes a land use area 311 with a land use field 312, a 
proportion controlled field 313, and land use scroll buttons 314. The land use 
scroll buttons allow the user to scroll through and add and delete land uses. The 
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land use field provides a drop-down list through which a user can select a land 
use of the selected subwatershed. The proportion controlled field 313 specifies 
the portion of the subwatershed that is controlled by the land use. A BMP field 
317 lists the BMPs that have been applied to the selected land use and provides 
a drop-down list of BMPs that can be applied. An initial year field 318 specifies 
the starting year of the associated BMP. A BMP area 316 lists the BMPs of the 
selected land use. The BMPs, shown in the BMP area at the same time, 
represent BMPs that are applied in series to the land use. BMP scroll buttons 319 
allow the user to scroll through sets of serial BMPs that are applied in parallel to 
the land use. The point source area includes a point source field 321, a treatment 
type field 322, and an installation year field 323. The point source field provides a 
drop-down list of the point sources of the selected subwatershed so that the user 
can select a point source to which a treatment is to be supplied. The treatment 
type field provides a drop-down list of the possible treatment types of the point 
source. The installation year field allows the user to specify when the selected 
treatment is installed. Point source buttons 324 allow a user to select, add, and 
delete different point sources. A view results button 331 allows a user to direct 
the WQTS to calculate the local credits of each subwatershed based on the 
specified controls (e.g., BMP and treatments). 
[0031] Figure 4 illustrates a report listing the local credits for a subwatershed in 

one embodiment. The report is for the currently selected subwatershed. The 
report lists future load with and without the controls applied, with their difference 
being the local credit. The report also summarizes the effect of the point source 
and non-point source controls. The summary of the non-point source controls 
includes an identification of the BMP, the install year, the life, the efficiency, the 
load reduction, the cost per pound of reduction, and the total cost of the control. 
The summary of the point source controls includes an identification of the point 
source, the treatment type, the effluent concentration, the install year, and the 
total cost of the control. This report can be generated at the subwatershed, 
watershed, and river basin levels depending on the selected calculation level. 
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[0032] After a user selects the river basin, the confluence point, the calculation 

level, and then the loading and controls buttons, the user selects the trading 
button to direct the WQTS to perform the trading calculations. The WQTS 
requests the user to select a confidence level and then calculates and displays 
the earned credits (aka the "watershed credits") at the selected confluence point 
for each subwatershed, factoring in the uncertainties and the selected confidence 
level. Figure 5A illustrates a graph of the earned credits for the subwatersheds 
contributing to a confluence point in one embodiment. The graph provides an 
indication of the credits versus year for each subwatershed. Although not 
illustrated in the figures the WQTS also allows a user to select pairs of 
subwatersheds for which trading ratios are to be calculated. Figure 5B illustrates 
a report of the earned credits for a subwatershed contributing to a confluence 
point in one embodiment. The report includes the local credits, the earned 
credits, and the uncertainty ratios for each year. 

[0033] Figure 6 illustrates a display page of a GIS interface for the WQTS in one 

embodiment. The display page 600 displays a map 601 of the subwatersheds 
contributing to a user-selected confluence point in a user-selected river basin. A 
layers area 602 allows a user to specify the information that is to be overlaid onto 
the map by selecting the appropriate boxes. For example, the user may select to 
overlay onto the map an outline of each subwatershed. After the trading 
information has been generated using the display page 100 of Figure 1, the 
WQTT toolbar 603 allows the user to display the results. Figure 7 illustrates a 
map overlaid with the loading information of the subwatersheds in one 
embodiment. Figure 8 illustrates a map overlaid with the local credit information 
of the subwatersheds in one embodiment. Figure 9 illustrates a map overlaid with 
the earned credit information of the subwatersheds in one embodiment. 

[0034] The following scenario provides an overview of a typical use of the WQTS. 

A person may need to offset an increase in the pollutant load of a point source. 
That person may have the option of adding certain BMPs to land uses under their 
control to offset the increase. Alternatively, it may be more cost-effective to apply 
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BMPs to land uses under the control of others. That person would use WQTS to 
enter the BMP information for their land use to generate the required earned 
credit to offset the increased load. That person would enter the BMP information 
and select a confidence level and be presented with cost and the earned credits 
with an associated uncertainty ratio. That person can then repeat the process for 
the BMPs that might be applied to others. That person can then use the trading 
ratios of the BMPs applied to their land uses and those of others to determine the 
equivalency of the earned credits. Assuming that equivalent credits can be 
obtained at a lower cost by applying BMPs to the land uses of others, then that 
person would attempt to trade credits with the other person who controls that land 
use. 

Loading Calculations 

[0035] The loading calculations performed by the WQTS in one embodiment are 

described in the following. 

[0036] The total current load of a pollutant for a subwatershed is calculated by 

summing the loads generated from each individual land use of the subwatershed. 
The load of a land use is calculated by multiplying the area of the land use in the 
subwatershed by the export coefficient of the land use. The export coefficient 
represents the amount of a pollutant generated by the land use in mass per unit 
area per time duration (e.g., pounds/acre/year). Thus, the total current load for a 
subwatershed is represented by the following equation: 

A=i(4*c,) (1) 

7=1 

where Li is the total current load of subwatershed /, Atj is the area of land use j in 
subwatershed /, Q is the export coefficient of land use 7, and ut is the number of 
land uses in subwatershed /. 
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[0037] 



The total future load of a pollutant for a subwatershed is calculated in a 
similar manner as for the total current load, except using future land uses. The 
total future load for a subwatershed is represented by the following equation: 



(2) 



where FLt is the total future load of subwatershed /, Aij is the area of future land 
use j in subwatershed /, Q is the export coefficient of land use j, and is the 
number of future land uses in subwatershed /. 
[0038] The standard deviation of the future load of a pollutant of a subwatershed 

is used to calculate future load at certain confidence levels. The standard 
deviation is the square root of the variance, which is calculated by summing the 
square of the product of future average load of each land use with the coefficient 
of variation for that land use. Thus, the standard deviation of the future load is 
represented by the following equation: 



1=1 



-ll/2 

(3) 



where o; is the standard deviation of the future load of subwatershed /, Aij is the 
area of future land use j within subwatershed /, Cj is the mean export coefficient 
of land use y, is the number of future land uses in subwatershed /, and CVCj is 
the coefficient of variability for land use j, 
[0039] The future load of a subwatershed at a certain confidence level is 

calculated by multiplying the future load by the standard deviation of the 
subwatershed by the z value of the confidence level. The future load at a certain 
confidence level is represented by the following equation: 

FLXx)^FLX\^z{x)cx,] (4) 
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where FLi(x) is the future load of subwatershed / at confidence level x, FU is the 
future load of subwatershed /, o; is the standard deviation of the future load of 
subwatershed /, and z(x) is the z value for confidence level x. 
[0040] The current point source load of a subwatershed is calculated by summing 

the loads of each point source in the subwatershed. The point source load is 
represented by the following equation: 

where PSU is the current point source load of subwatershed /, Z,y is the current 
load of current point source j in subwatershed /, and /w, is the number of current 
point sources in subwatershed /. 
[0041] The future point source load of a subwatershed is calculated in a similar 

manner to the current point source load and is represented by the following 
equation: 

FPSL,=^L, (6) 

where FPSLi is the future point source load of subwatershed /, Ly is the future load 
of future point source j in subwatershed /, and /w, is the number of future point 
sources in subwatershed /. 
[0042] The delivery ratio represents the ratio of the load of a subwatershed that 

makes it downstream to a confluence point. Each subwatershed has a local 
delivery ratio that represents the portion of the inflow load that makes it to the 
outlets of the individual watersheds. The delivery ratio is calculated by multiplying 
the local delivery ratio of each downstream subwatershed together and is 
represented by the following equation: 

DR,=Y\LDR^ (7) 
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where DRi is the delivery ratio of subwatershed /, LDRj is the local delivery ratio of 
downstream subwatershed y, and w, is the number of downstream subwatersheds 
of subwatershed 7, including subwatershed / itself. 
[0043] The point source load of a subwatershed that makes it to the confluence 

point is the point source load of the subwatershed multiplied by the delivery ratio. 
The point source load of a subwatershed contributing to a confluence point is 
represented by the following equation: 

PSi^PSL^^DRi (8) 

where PSi is the point source load contribution of subwatershed / at the 
confluence point, where PSLt is the current point source load of subwatershed /, 
and DRi is the delivery ratio of subwatershed /. 
[0044] The future point source load of a subwatershed that makes it to the 

confluence point is calculated in a manner similar to the way the current point 
source load is calculated and is represented by the following equation: 

FPS,=FPSL,^DR, (9) 

where FPSt is the future point source load contribution of subwatershed / at the 
confluence point, FPSLt is the future point source load of subwatershed /, and DRi 
is the delivery ratio of subwatershed / to the confluence point. 

Local and Earned Credit Calculations 
[0045] The local and earned credit calculations performed by the WQTS in one 

embodiment are described in the following. The non-point source load from a 
subwatershed having at most one BMP per land use is calculated based on the 
removal efficiency of the BMP and the portion of the land use affected by the BMP 
and is represented by the following equation: 
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[0046] 



where Ly is the load of subwatershed / from land use Q is the export coefficient 
of land use j, Ay is the area of land use j in subwatershed /, and ^ is the fraction of 
land use j in subwatershed / treated with a BMP having a pollutant removal 
efficiency of ry. 

The mean (or average) load per interval (e.g., season) of a subwatershed 
over a period (e.g., one year) is calculated based on the mean of the export 
coefficient and the removal efficiency. The mean load of a subwatershed is 
represented by the following equation: 



= c, * 4 * (1 - ^/ ) * 4 + * 4 * (1 - 4. ) 



(11) 



[0047] 



where Z,y is the mean load of subwatershed / from land use j, Cj is the mean 
export coefficient of land use j, Ay is the area of land use j in subwatershed /, and 
fy is the fraction of land use j in subwatershed / treated with a BIVIP having a mean 
pollutant removal efficiency of fjj . 

The coefficient of variation of the load of a subwatershed is derived based 
on the coefficients of variation of the export efficiency and the removal efficiency 
and is represented by the following equation: 



CVL,= 



(l + CVC/) 



1/2 



(12) 



where CVLy is the coefficient of variation of the load generated from land use j 
within subwatershed 1, CVCj is the coefficient of variation of the export coefficient 
of land use 7, CVrjj is the coefficient of variation of the removal efficiency of the 
BMP, and /y is the fraction of land use j within subwatershed / treated with the 
BMP having a mean pollutant removal efficiency of rfj . 
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[0048] Equations (11) and (12) for the mean and coefficient of variation of the land 

use assume implementation of only one BMP per land use as a pollution 
controlling unit. When multiple BMPs are used in series on a land use, then the 
removal efficiency is calculated as one minus each of the inefficiencies (i.e., I-tj) 
multiplied together. The mean removal efficiency of the group of multiple BMPs in 
series combination is represented by the following equation: 

^=i-n(i-^o <^3A) 

where fj^ is the mean removal efficiency of the serial BMPs of land use j in 

subwatershed /, 7,^ is the mean removal efficiency of BMP k in land use j, and «/, 

is the number of BMPs of land use j in subwatershed /. Using equations (13A) 
and (11), the mean pollutant loading from a land use having a group of multiple 
BMPs in series is represented by the following equation: 

where Lj is the mean load of subwatershed / from land use 7, Q is the export 
coefficient of land use 7, Ay is the area of land use j in subwatershed z, fy is the 
fraction of land use j in subwatershed / treated with a BMP having a mean 
pollutant removal efficiency of ifj , and riy is the number of BMPs of land use j in 

subwatershed /. 

[0049] The coefficient of variation of the removal efficiency for serial BMPs is the 

product of one plus the variance of the removal efficiency divided by the 
inefficiency for each BMP. That coefficient of variation is represented by the 
following equation: 
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AI2 



(14A) 



where CF% is the coefficient of variation of the removal efficiency of land use j in 
subwatershed /, Var(Tjjif) is the variance of the removal efficiency for BMP k in land 
use j, fjji is the mean removal efficiency for BMP k in land use j, and nu is the 

number of BMPs of land use J in subwatershed /. Using equations (14A) and (12), 
the coefficients of variation of pollutant loading from land use j having a group of 
multiple BMPs in series are represented by the following equation: 



CVL^.= 



(i+cFc;) 



fim-^^) n 

k=\ k=\ 



1 + - 



i-4+^n(i-^*) 



-1 



1/2 



(14B) 



where CVLij is the cx)efficient of variation of the load of subwatershed /, CVQ is the 
coefficient of variation of the export coefficient of land use 7, Ay is the area of land 
use j in subwatershed z, ffj^ is the mean removal efficiency of BMP k in land use 7, 

and riij is the number of BMPs of land use j in subwatershed /. 
[0050] When multiple BMPs are in parallel on a land use, then the mean load for a 

land use in the subwatershed is the product of the mean export coefficient of that 
land use for the subwatershed times the sum of the areas of the land use adjusted 
by the mean removal inefficiency. The load with multiple BMPs in parallel is 
represented by the following equation: 
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where L^j is the mean load of subwatershed / from land use y, Cj is the mean 

export coefficient of land use y, Aij is the area of land use j in subwatershed /, //^^ 

is the mean removal efficiency for BMP k in land use 7, and riij is the number of 
BMPs in parallel of land use j in subwatershed /. 
[0051] The coefficient of variation of pollutant loading for BMPs in parallel is 

represented by the following equation: 



.1/2 



(l+CFC/) 



1 — 



-1 



J J 



(16) 



where CVLij is the coefficient of variation of the load of subwatershed / from land 
use 7, CVCj is the coefficient of variation of the export coefficient of land use 7, Ay is 
the area of land use j in subwatershed /, ff^^ is the mean removal efficiency of 

BMP k in land use y, CVrjjk is the coefficient of variation of BMP k in land use y, and 
«y is the number of BMPs of land use y in subwatershed /. 
[0052] The mean pollutant load of a subwatershed is the sum of the mean 

pollutant loads from each land use within the subwatershed and is represented by 
the following equation: 

A=sr, (17) 

where Z, is the mean pollutant load of subwatershed /, Lij is the mean pollutant 
load of land use j in subwatershed /, and w, is the number of land uses in 
subwatershed /. 
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[0053] The variance of the total load from a subwatershed is represented by the 

following equation: 



Varm) = f^(L^*CVLyf (18) 



where Var(Li) is the variance of the pollutant load of subwatershed /, L. is the 

mean pollutant load of subwatershed z, CVLij is the coefficient of variation of the 
load of land use j in subwatershed z, and is the number of land uses in 
subwatershed z. 

[0054] The coefficient of variation of the total load from a subwatershed is 

represented by the following equation: 



\\/2 



^ (MAIL (19) 

where CVLt is the coefficient of variation of the load of subwatershed /, VarfLi) is 
the variance of the pollutant load of subwatershed /, and is the mean load of 
subwatershed /. 

[0055] The one-sided lower and upper confidence intervals for an assumed 

distribution of the load delivered from a subwatershed at its outlet is represented 
by the following equation: 

a,=Z;±z,*Z;*[Far(l,)f =A[i±^.*^ (20) 

where CI. is the confidence interval for the load delivered from the outlet of 

subwatershed z, L. is the mean load of subwatershed z, Var(Li) is the variance of 

the pollutant load of subwatershed z, and is the z value at significance level a . 

[0056] The two-sided lower and upper confidence interval for an assumed 

distribution of the load delivered from a subwatershed at its outlet is represented 
by the following equation: 
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c/,=A.±z„/3* 



Z;.*[Far(A.)r=A{l±-„.*[^«KA)r] 



(21) 



where C/,. is the confidence interval for the load delivered from the outlet of 
subwatershed /, ^ is the mean load of subwatershed i, Var(L^ is the variance of 
the pollutant load of subwatershed /, and z^/^ the z value at significance level 
all. 

The uncertainty ratio for a non-point source accounts for the variability in 
load delivered from a subwatershed due to variation in the export coefficient of a 
land use and the removal efficiency of a BMP. The uncertainty ratio is calculated 
by dividing the mean load by the mean load times one plus the z value at a 
certain confidence level (e.g., 95%) times the coefficient of variation of the load. 
When the value of Za corresponding to the 95*" percentile of the total load 
distribution is approximately 2, the uncertainty ratio for the non-point source load 
delivered from a subwatershed is represented by the following equation: 



where URi is the uncertainty ratio for subwatershed /, L. is the mean load of 
subwatershed /, CVLt is the coefficient of variation of the load of subwatershed /, 
and Za is the z value at a certain confidence level. 

The uncertainty ratio for a point source, such as a wastewater treatment 
plant, is assumed to have negligible variabiHty in the load delivered to the outlet of 
a subwatershed (i.e., CV \s equal to 0). Therefore, the uncertainty ratio due to the 
point sources of a subwatershed Is represented by the following equation: 



where C/R, is the uncertainty ratio of a point source in subwatershed /. 

The exchange ratio for a given point or non-point source in a watershed / is 
represented by the following equation: 



UR,= 




(22) 



L\\ + z^*CVL,\ \ + 2*CVL, 



UR,=\ 



(23) 
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ER, = DR, * UR, * WR, * RR, (24) 

where ERi is the exchange ratio for subwatershed /, DRi is the delivery ratio for 
subwatershed /, URi is the uncertainty ratio for subwatershed /, WRi is the water 
quality ratio for subwatershed /, and RRi is the retirement ratio for subwatershed /. 
The water quality ratio is included to protect local water quality from adverse 
impacts of pollutant loading contributed by the subwatershed as it travels to the 
confluence point. The retirement ratio allows a certain percentage of the 
available credits to be retired or excluded from trading to increase the potential for 
water quality benefit or to provide a margin of environmental safety for the overall 
trading program. 

[0060] Since the delivery ratio as calculated above does not take into 

consideration any variability, a delivery ratio correction factor may be used to 
adjust the delivery ratio. The correction factor is assumed to have a normal 
distribution with a unity mean and a certain coefficient of variation based on 
uncertainty in the delivery ratio. The adjusted exchange ratio is represented by 
the following equation: 

ERi=X^ER, (25) 

where ERi is the corrected exchange ratio for subwatershed / and X is the 
correction factor. 

[0061] The trading ratios are represented by the following equation: 

FR 

TR.-^ (26) 

where TRs is the trading ratio for seller s, ERs is the exchange ratio for seller 5, and 
ERb is the exchange ratio for buyer b, 
[0062] The amount of credits earned by a subwatershed by applying a BMP is 

represented by the following equation: 
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ECi=(Lr-Lr')*ERi 



(27) 



where Ed is the mean earned credit for subwatershed i, Lj"" is the mean load of 
subwatershed / prior to applying the BMP, Lj^ is the mean load of subwatershed 

/ after applying the BMP, and ERi is the corrected exchange ratio for 
subwatershed /. 

[0063] The variance of the earned credit is represented by the following equation: 

Var (jECi ) = ((Zf" * CVLf" )' + [Zr' * CVLf^ )' ) * (28) 

where Far(£C, ) is the variance of mean earned credit for subwatershed /, Lj'" is 

the mean load of subwatershed / prior to applying the BMP, CVir is the 
coefficient of variation of subwatershed / phor to applying the BMP, Lf°" is the 
mean load of subwatershed i after applying the BMP, CVir" is the coefficient of 
variation of subwatershed / after applying the BMP, and ERi is the corrected 
exchange ratio for subwatershed /. 
[0064] The lower and upper confidence intervals of the earned credit 

conresponding to a certain significance level are represented by the following 
equation: 

C/, =£C,[l±z.„[Far(£C,)]"'] (29) 



where C/,. is the confidence interval of subwatershed /, EQ is the mean earned 

credit for subwatershed /, Za/2 is the z value at significance level a/2, and VarfEQ 
is the earned credit for subwatershed /. 
[0065] Figures 10-14 illustrate the processing of the WQTS in one embodiment. 

Figure 10 is a block diagram illustrating the components of the WQTS in one 
embodiment. The WQTS includes a user interface component 1001, a GIS 
component 1002, a loading component 1003, a controls component 1004, a 
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trading component 1005, and a database 1006. The user interface component 
controls the user interface of the WQTS. The GIS component provides a 
geographic information system interface of the WQTS. The loading component 
performs the loading calculations of the WQTS including the current loads and 
future loads. The controls component controls the input of the control information 
and performs the controls calculations of the WQTS including the local credits. 
The trading component performs the trading calculations of the WQTS including 
the earned credits. The database contains the GIS information and the other 
information describing BMPs, coefficients of variation, and so on. 
[0066] Figure 1 1 is a flow diagram illustrating the processing of a component that 

identifies subwatersheds in one embodiment. After a river basin and a confluence 
point are selected, the component identifies the subwatersheds that are upstream 
from the confluence point. The illustrated processing is recursive. One skilled in 
the art will appreciate that non-recursive processing may alternatively be 
performed. The component is passed an identification of the confluence point. In 
blocks 1 101-1 104, the component loops selecting each subwatershed that directly 
contributes to the confluence point and recursively identifies each subwatershed 
that indirectly contributes to the confluence point. In block 1101, the component 
selects the next subwatershed that directly contributes to the passed confluence 
point. In decision block 1102, if all the subwatersheds have already been 
selected, then the component returns, else the component continues at block 

1103. In block 1103, the component adds the selected subwatershed to a list of 
subwatersheds that are upstream from the original confluence point. In block 

1 104, the component recursively invokes the component passing the inflow of the 
selected subwatershed as a confluence point to identify the subwatersheds that 
contribute to the selected subwatershed. The component then loops to block 
1101 to select the next subwatershed that directly contributes to the passed 
confluence point. 

[0067] Figure 12 is a flow diagram illustrating the processing of a component that 

calculates the loading for the subwatersheds that are upstream from a confluence 
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point in one embodiment. The component calculates the current and future loads 
for the point and non-point sources and calculates the future loads of the non- 
point sources at various confidence levels. In block 1201 , the component selects 
the next subwatershed. In decision block 1202, if all the subwatersheds have 
already been selected, then the component completes, else the component 
continues at block 1203. In block 1203, the component calculates the current and 
future loads for the non-point sources of the selected subwatershed as illustrated 
by equations (1) and (2). In block 1204, the component retrieves the future load 
of the selected subwatershed from the database. In block 1205, the component 
calculates the future load of the non-point sources of the selected subwatershed 
at various confidence levels as illustrated by equations (3) and (4). The 
component then loops to block 1201 to select the next subwatershed. 
[0068] Figure 13 is a flow diagram posted in the processing of a component that 

calculates the local credits for the subwatersheds that are upstream from a 
confluence point in one embodiment. The component loops selecting each 
subwatershed and calculating the local credits for point sources and non-point 
sources of that subwatershed. In block 1301, the component selects the next 
subwatershed. In decision block 1302, if all the subwatersheds have already 
been selected, then the component returns, else the component continues at 
block 1303. In block 1303, the component calculates the load of the point 
sources contributing to the local credits of the selected subwatershed. In decision 
blocks 1304-1306, the component loops calculating the load of each land use 
contributing to the local credits of the selected subwatershed. In block 1304, the 
component selects the next land use of the selected subwatershed. In decision 
block 1305, if all the land uses of the selected subwatershed have already been 
selected, then the component continues at block 1301 to select the next 
subwatershed, else the component continues at block 1306. In block 1306, the 
component evaluates the BMPs of the selected land use. The component 
calculates the reduction in load attributable to the serial BMPs according to 
equations (10) and (12) and the parallel BMPs according to equations (10), (14A), 
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and (14B) to identify the portion of the load contributed by that land use. The 
component then loops to block 1304 to select the next land use. The component 
adds the load for each land use to determine a total non-point source load. The 
component then adds the load for the point sources and the non-point sources 
and subtracts the sum from the future load calculated without the controls to 
determine the local credits for the selected subwatershed. 

[0069] Figure 14 is a flow diagram illustrating the processing of a component that 

calculates the earned credits for subwatersheds that are upstream from a 
confluence point in one embodiment. In block 1401, the component selects the 
next subwatershed. In decision block 1402, if all the subwatersheds have already 
been selected, then the component returns, else the component continues at 
block 1403. In block 1403, the component calculates the coefficient of variation of 
the load for the selected subwatershed according to equation (18). In block 1404, 
the component calculates the uncertainty ratio for the selected subwatershed 
according to equation (22). In block 1405, the component calculates the 
exchange ratio for the selected subwatershed according to equations (24) and 
(25). In block 1406, the component calculates the earned credits for the selected 
subwatershed according to equation (27). The component then loops to block 
1401 to select the next subwatershed. 

[0070] From the foregoing, it will be appreciated that specific embodiments of the 

invention have been described herein for purposes of illustration, but that various 
modifications may be made without deviating from the spirit and scope of the 
invention. One skilled in the art will appreciate that the WQTS can be used to 
evaluate for various types of pollutants (e.g., phosphorus and sediment) and can 
be used at various levels of granularity (e.g., subwatersheds, watersheds, and 
point and non-point sources). Accordingly, the invention is not limited except as 
by the appended claims. 
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